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Abstract 

The long range structure of the nucleon is discussed 
starting from the old model of a quark bag with a pion 
cloud ("cloudy bag") carrying on to the more recent 
ideas of the parton model of the nucleon. On the 
basis of the most recent measurements of the form 
factors at MAMI, JLab and MIT quantitative results 
for nucleon charge densities are presented within both 
non-relativistic and relativistic frameworks. 

1 Introduction: Ranges 

Many physicists have the following picture of the 
structure of the nucleon : In the inner region at "short 
range" reside quarks bound by gluons, in the outer 
region at "large distances" live mesons and in partic- 
ular pions. The nucleon consists of constituents, i.e. 
"constituent quarks" and pions, as the atom consists 
of a nucleus and electrons, and the nucleus of protons 
and neutrons. The range of the inner "coloured" re- 
gion, frequently called "confinement radius", is rather 
elusive. It is a model parameter in the old bag model 
or models for the nucleon resonances based on the 
constituent quarks. One can estimate it from ex- 
periment by identifying it with the "annihilation ra- 
dius" of the antiproton-proton system. Only quark- 
antiquark annihilation in the overlap region of colour 
can contribute to annihilation [lj. It amounts to ap- 
proximately 0.8 fm, in reasonable agreement with the 
mentioned model parameters. However, it cannot be 
easily identified with the root-mean-square (rms) ra- 
dius of the electric charge of the proton, since the 
pion cloud will contribute to the charge distribution. 
A rough idea of the range of this contribution may be 
gotten by the Compton wave length of the pion which 
is of the order of 1.4 fm. A similar picture emerges 



from diffractive scattering of high energy protons and 
antiprotons [2J. 

However, this picture fails in two ways. Firstly, 
the nucleon moves after the scattering in most ex- 
periments at relativistic velocities and therefore its 
structure looks different in different reference frames. 
The simplest example is the transformation of the 
magnetic moment into an electric dipole moment. 
This situation makes it particularly difficult to com- 
pare experiments to model calculations based on rest 
frame wave functions. These wave functions have to 
be "boosted" to the correct momentum transfer and 
there is no consistent way of doing that. Secondly, at 
relativistic energies particle-antiparticle, i.e. quark— 
antiquark, pairs have unavoidably to be considered 
making the picture much more involved. These two 
aspects will be discussed in Section[3] 

The relativity destroys the simple picture also in an- 
other way. We usually relate ranges with momentum 
transfer via Heisenberg's uncertainty relation. How- 
ever, since in relativistic mechanics space and time are 
intimately connected no relativistic uncertainty rela- 
tion exists [3J. Therefore the assignment of ranges 
to quantities depending on the negative four momen- 
tum transfer squared Q 2 = —q 2 , as e.g. in some 
plots of the running coupling constant, is rather mis- 
leading. We shall, therefore, in this article distinguish 
between the non-relativistic picture derived at small 
momentum transfers and the relativistic case where 
we have to use the relativistic quantum field theo- 
retic description. At small momentum transfers we 
can approximate the four momentum transfer q with 
the three momentum transfer q 2 rj — q 2 and maintain 
the familiar interpretation of form factors (Section^. 
In Section[3] we shall show how we can connect the 
non-relativistic picture to the underlying quark-gluon 
structure. We shall see that new experiments in just 
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the relativistic domain are needed in order to clarify 
how nucleons are made up of quarks and gluons or 
more precisely how hadrons emerge from QCD. 

2 Nonrelativistic interpretation of 
nucleon form factors 

The non-relativistic electromagnetic form factor (FF) 
is a special case of one of the fundamental observables 
in quantum mechanics. It is the matrix element of the 
interaction propagator of the exchange photon, i.e. 
the operator e lq ' r . It appears in all domains of micro 
physics ranging from atomic physics, the Mossbauer 
effect to particle physics |4j. It reads in the case of 
elastic scattering from a charge density p{f) = (f)r) 



r (r e 



iq-r 



\r). 



(1) 



Since we have the idea that the photon interacts with 
a single constituent we can interpret the form factor 
as the probability that we scatter from a constituent 
with a momentum just right to leave the bound sys- 
tem "intact", i.e. in the ground state. Since one 
deals mostly with spherically symmetric systems it is 
customary to write F(q 2 ). 

The method to derive charge and magnetic distri- 
butions from the scattering of electro-magnetically in- 
teracting particles - almost exclusively these are elec- 
trons - has been extensively used for nuclei and nu- 
cleons. As a reminder we present the Rosenbluth for- 
mula, the key formula connecting the measured cross 
sections to the FFs for a spin-1/2 particle : 
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where E is the incoming and E' the outgoing electron 
energy, 9 the scattering angle, and r = Q 2 /(AM 2 ) 
with M the mass of the nucleon. Dividing the mea- 
sured cross section by the Mott cross section, i.e. the 
cross section for point particles, yields the combina- 
tion of the quadratic electric form factor G\ and mag- 
netic G 2 M . By varying the scattering angle 9 and the 
incoming electron energy E in such a way that Q 2 



is constant, G 2 E and G 2 M can be separated, the so 
called Rosenbluth separation. This formula is rela- 
tivistically correct meaning that the FFs are functions 
of Q 2 . As a practical tool, the Rosenbluth method is 
not ideal, however, for very small cross sections as is 
the case when extracting the small G 2 E term of the 
neutron, or when separating at large Q 2 the small G 2 E 
contribution of the proton from that of G\ [ since the 
latter dominates at large momentum transfers. Here 
polarized electrons available with sufficient intensity 
since about a decade, together with polarized targets 
or recoil polarimetry, have changed our possibilities 
dramatically. 

In order to give a feeling for the significance of the 
method we present the example of nuclei. Nuclei are 
heavy and have in many cases no magnetic contribu- 
tion, hence the charge distribution can be precisely 
determined. Figure[T] shows a collection of data for 
the scattering of electrons from le O and 20S Pb. 

The Fourier transform (see Eq. ([!])) of the FFs de- 
rived from these cross sections allows for a model in- 
dependent precise determination of the charge dis- 
tributions. This is shown in Fig. [2] for some nuclei. 
These distributions can be well modeled by mean field 
calculations based on the nuclear shell model [6J. We 
see that the nuclear charge distributions have a subtle 
structure on top on the bulk behaviour of the satu- 
rated nuclear matter. This structure reflects the shell 
model wave functions and is very sensitive to the de- 
tails of the nucleon-nucleon interaction but even more 
importantly to the many body features of the nuclei 
made up of the constituents protons and neutrons. 
Only the most sophisticated mean field theories are 
able to almost reproduce this structure [6]. It re- 
mains an unexplained deviation which may be due to 
too narrow error bands [7j. 

Let us now turn to the nucleons and see what we 
know here (see e.g. Refs. [8] [9| [10] for some recent 
reviews on nucleon form factors). The pre-2000 data 
suggest that the magnetic and electric form factor of 
the proton follow a universal form, the "dipole form" 
G D (Q 2 ) = l/(l+Q 2 /A 2 D ) 2 , with the scale parameter 
A E = 0.843 GeV approximately equal to the mass m p 
of the p meson. The case was closed and considered 
to be text book material. This finding was the ba- 
sis for the much discussed "vector dominance model" 
in strong interactions. (For a recent discussion see 
ref. The p meson mass defines a "small range" 

of an exponential distribution, somewhat unphysical 
due to its discontinuity at r = 0. It was believed that 
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Figure 1: A compilation of form factors (FFs) from 
elastic electron scattering for le O and 208 P6 mea- 
sured at the electron-accelerator laboratories: Am- 
sterdam, Darmstadt, Mainz, Saclay, and Stanford. 
The FFs have been corrected for the Coulomb dis- 
tortions different at the electron energies available 
at these laboratories. The diffraction pattern deter- 
mining the charge distribution is nicely seen. (From 
Ref. 0). 



the dipole form would also describe the long ranges 
characterized by the rms radius of the proton. A 
measurement at the High Energy Physics Laboratory 
HEPL at Stanford gave y/ir^j = (0.805 ± 0.011) fm 
~ \jYljm p |12| . However, as it is now clear, a 
measurement at Mainz in 1980 [13J which yielded 
yfJrZ) = (0.862 ± 0.012) fm was closer to the best 
value (0.882 ± 0.010) fm available today. (For a more 
detailed discussion see Ref. |14j.) This value is con- 
firmed by a recent high precision determination at the 
Mainz Microtron (MAMI) yielding [TJ: 

ypy= (0.879 ±0.005 sta t.± 

0.004 syst . ± 0.005 modcl ) fm. (3) 

As discussed in Ref. [14] this value is significantly 




— Mean-Field 

Theory 
Experiment 



Figure 2: Charge distributions of spherical nuclei as 
derived by model independent analyses from exper- 
iment (solid lines). The hatched areas indicate the 
range of the charge consistent with the data. The 
dashed lines depict the results of calculations in 
the framework of nuclear mean-field theory. (From 
Ref. [6]) 



larger than the largest value derived from dispersion 
relations and is unexplained in most nucleon mod- 
els. On the other hand, a recent study of the Lamb 
shift in muonic hydrogen at the Paul-Scherrer Insti- 
tute in Zuerich yielded a very precise value as low as 
(0.840 ± 0.001) fm [15J in agreement with the upper 
bound of the dispersion-relations calculations. It is 
not credible to assign this five standard deviation dif- 
ference to a deficiency in the theory of the electronic 
experiments which is the same Quantum Electrody- 
namics as for the muonic hydrogen. However, the 
calculations of the Lamb shift in muonic hydrogen 
are difficult due to the strong distortions of the muon 
wave function calling for further theoretical work. 

Considering the focus on "long ranges" the ques- 
tion arises how much of the rms radius is possibly due 
to a pion cloud. There were weak indications that 
the pion cloud could be directly seen in the FFs in a 
similar way as the shell structure of the nucleus. This 
indication came from a coherent analysis of the data 
available until 2003 for the electric and magnetic FFs 
of the proton and neutron |16|. Since 2003 the data 
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base has improved so much that we want to base this 
discussion on the most recent measurements. 

Figure[3] shows the electric FF of the proton as de- 
rived from a direct fit of a FF model to data ob- 
tained with the 3-Spectrometer set-up at MAMI. This 
method takes advantage of modern computers and 
fits Eq.Q to a large set of angular distributions 
measured at five energies 180, 370, 450, 720, and 
850 MeV. All together about 1400 settings were mea- 
sured. In this way the "measurement at constant Q 2 " , 
i.e. the old Rosenbluth separation, becomes obsolete 
and a very broad kinematic range can be covered in- 
deed. However, since this method is somewhat unfa- 
miliar Fig.|4]shows the best fit curve of Fig.[3]together 
with the results at those kinematics at which a tra- 
ditional Rosenbluth separation could be done. The 
agreement is very good. A point of concern may be 
the analytical model used for the electric form factor 
Ge p and magnetic form factor Gmp in the direct FF 
fits. Here about a dozen different forms have been 
used all yielding essentially the same results [7j. The 
rms radius is, however, somewhat dependent on them 
and the second systematic error in Eq. ^ reflects this 
dependence. 

It is evident that these form factors show some 
structure after the gross dependence, assumed to be 
given by the standard dipole form, has been divided 
out. In Fig. [3] one observes two slopes for Ge/Gd- 
The steep negative slope at small Q 2 is reflected in 
the large rms radius discussed above. The reverse 
is true in Fig. [5] for the rms radius of Gu /(Mp Gd)- 
A shoulder structure is indicated in both FFs at 
Q 2 s» 0.15 GeV 2 . It is shifted compared to the bump 
structure derived by Friedrich and Walcher [16J from 
the pre-2003 data and cannot be identified with it. 
However, just considering the scale of the rms radius 
and the scale of the structure it is suggestive to look 
for pion cloud contributions in the modelling of the 
nucleon. 

Since the bulk charge of the proton resides at "small 
ranges" and extends out to the range of the pion 
cloud, the separation of the inner component from 
the pion cloud contribution will be somewhat arbi- 
trary. Here the neutron with a total zero bulk charge 
promises an experimental access since one way the 
neutron could acquire a charge distribution is just by 
its virtual dissociation n — > p + ir~ . This means that 
the pion cloud should be more clearly visible again as 
a signal at large radii r w Xcompton = « 1.4 fm 

in the neutron charge distribution. 
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Figure 3: The electric FF of the proton Ge p /Gd 
obtained from a direct fit to the cross sections 
with the spline model to data measured with the 
3-Spectrometer set-up at MAMI. The FF is nor- 
malised to the dipole form given in the text. The 
1 (T-error band is shown for the indicated errors. In 
order to show the improvement, the old data ob- 
tained with the classical Rosenbluth separation are 
also depicted. 



One could be tempted to look at the neutron rms 
radius as in the case of the proton in first place. How- 
ever, the mentioned recoil effect causes the magnetic 
moment to contribute to the small electric form fac- 
tor. It turns out that the electric rms radius of the 
neutron is a subtle interplay between the recoil ef- 
fect and the charge radius proper. We do not want 
to elaborate this here. A summary can be found in 
Ref. [14j. Instead we will show the most recent results 
for the electric FF of the neutron. 

Figure|6] shows a compilation of all data including 
the results from the MIT Bates measurement with 
the BLAST detector at the South Hall Ring at small 
Q 2 [18J and the measurement at the Jefferson Lab- 
oratory (JLab) by the Hall A Collaboration at large 
Q 2 [19]. All these measurements use the depen- 
dence of the cross sections of polarized targets or 
recoil polarizations of the ejected nucleons for polar- 
ized electrons from the FFs. The curve shows a fit of 
the phenomenological model of Friedrich and Walcher 

(fw) ng. 

As already mentioned, the phenomenological FW 
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Figure 4: As Fig. [3] without error bands but with 
the results from the classical Rosenbluth separation 
of the new MAMI data. 

model tried a coherent fit of a smooth bulk curve with 
a superimposed bump for the pion cloud. This fit of 
the before 2003 data showed indeed a 2 a contribution 
above the smooth curve causing a shift of the nega- 
tive charge to radii around 1.5 fm. The inclusion of 
the MIT [18] and JLab data [19] in the fit makes the 
bump disappear. At small Q 2 the neutron FFs cannot 
be reliably determined due to the model dependence 
of the extraction of the form factor from the scatter- 
ing of the polarized electrons from unavoidably bound 
neutrons in deuterium or 3 He targets. Only new even 
more precise measurements will be able to improve 
the situation. However, it will be mandatory to fur- 
ther study the reaction mechanism experimentally in 
order to check the theoretical corrections. 

The scope of this article is on "long range struc- 
ture" meaning that one wants to discuss the idea of 
the spatial distribution of the constituents of the nu- 
cleon. In the rest frame the three-dimensional charge 
distribution of a spherically symmetric non-relativistic 
system is obtained as : 

p 3d (r) = J ^k 2 j (rk)p(k), (4) 

where p(k) is an intrinsic FF. As already pointed out 
the relativistic effects do not allow for a simple inter- 
pretation of the electric and magnetic FFs in terms 
of charge and magnetic density distributions in the 



Figure 5: As Fig. [3] for the magnetic FF 
GMp/il^pGn) of the proton. 

rest frame system. However, there is special refer- 
ence system, the Breit or brick-wall system, defined by 
having no energy transfer v to the nucleon, in which 
the charge operator for a non-relativistic (static) sys- 
tem is only expressed through the electric FF Ge- 
In this system the four-momentum transfer squared 
becomes q 2 = v 2 — q 2 = —q 2 . However, the rest 
frame systems (laboratory systems) for different q 2 
move with different velocities with respect to the Breit 
system. For a relativistic reference system, to relate 
its intrinsic FF p(k) and density to the Breit frame 
in which the system of mass M moves with veloc- 
ity v = WrJJl + t), requires a Lorentz boost relat- 
ing k 2 = Q 2 /(l +t). This relation shows that for 
Q 2 — > oo, there is a limiting largest intrinsic wave 
vector k —> 2M = 2w/Xu m . In the rest frame, no in- 
formation can be obtained on distance scales smaller 
than this wavelength due to relativistic position fluc- 
tuations (known as the Zitterbewegung) . For a non- 
relativistic system as le O, displayed in Fig. [5] Xu m 
is below 0.04 fm, whereas for the nucleon, it corre- 
sponds with Xum ~ 0.66 fm. Extracting the density 
for a relativistic system as the nucleon, therefore re- 
quires a prescription in order to relate the intrinsic 
FFs p(k) in Eq. Q to the experimentally measured 
FFs. As an example, in Ref. [20) . the prescription 
p(k) = (1 + T) n G E (Q 2 ) motivated by simple mod- 
els of nucleon structure was used and parameter val- 
ues n = 0,1, and 2 were investigated. To see the 
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Figure 6: A compilation of the recent data for the 
neutron electric FF, Geu, obtained at NIKHEF 
(triangle), Bates MIT (stars), JLab (squares), and 
MAMI Mainz (circles), with the quasi free ed —> 

pfte and the e 3 He — > ppne reactions. Either the 
polarization of th e recoiling neutron ft is measured 

or the target 3 He, and with it the neutron, is polar- 
ized. The protons p and pp are unobserved, respec- 
tively. The curves correspond with a new fit of a 
phenomenological model (solid red curve) [TB] , and 
of a Generalized Parton Distribution parametriza- 
tion (dotted black curve) [T7] . 



transitional region from the distance scales where rel- 
ativists position fluctuations hamper our extraction 
of rest frame densities to distances where the con- 
cepts of a non-relativistic many-body system can be 
approximately applied, we visualize the charge density 
in Fig. [7] It depicts the Fourier transform Eq. Q us- 
ing the fit in Fig. [6] (solid red curve). One notices a 
negative charge density at distances around and larger 
than lfm. With all caveats we may interpret the neg- 
ative charge as a "pion cloud" in the nonrelativistic 
limit since it extends beyond the confinement radius 
of about 0.8 fm. 



3 Relativistic picture 

We now turn to the relativistic picture and see how 
it does complicate matters, however, for the benefit 
of a deeper insight. As already mentioned both, the 
size and the shape of an object, are not relativistically 
invariant quantities: observers in different frames will 
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Figure 7: Charge distribution of the neutron as 
derived from the Fourier transform of the Geu fit 
(solid red curve in Fig. [6]) . The dashed part of the 
curve is for r < Xu m — 2it/(2M), where one is in- 
trinsically limited to resolve the density due to the 
Zitterbewegung of the nucleon. 



infer different magnitudes for these quantities. Fur- 
thermore when special relativity is written in a co- 
variant formulation, the density appears as the time 
component (zero component) of a four-current den- 
sity J M = (p, J) (in units in which the speed of light 
c=l). 

Besides the relativistic kinematic effects, as e.g. the 
length contraction, the concept of size and shape in 
relativistic quantum systems, such as hadrons, is also 
profoundly modified as the number of degrees of free- 
dom is not fixed anymore. In relativistic quantum 
mechanics the number of constituents of a system is 
not constant as a result of virtual pair production. 
We consider as an example a hadron such as the pro- 
ton which is probed by a space-like virtual photon, 
as shown in Fig.[8j A relativistic bound state as the 
proton is made up of almost massless quarks. Its 
three valence quarks making up for the proton quan- 
tum numbers, constitute only a few percent of the 
total proton mass. In such a system, the wave func- 
tion contains, besides the three valence quark Fock 
component \qqq), also components where additional 
qq pairs, the so-called sea-quarks, and (transverse) 
gluons g contribute leading to an infinite tower of 
\qqqqq), \qqqg) , ... components. When probing such 
a system using electron scattering, the exchanged vir- 
tual photon will couple to both kind of quarks, valence 
and sea, as shown in Figs. [8] (a) and (b). In addition, 
the virtual photon, can also split into a qq pair, lead- 
ing to a transition from a 3q state in the initial wave 



G 




Figure 8: Coupling of a space-like photon to a rela- 
tivistic many-body system, such as the proton. Top 
panel (a) : diagonal transition where the photon 
couples to a quark in the leading 3q Fock compo- 
nent of the proton. Middle panel (b) : diagonal 
transition where the photon couples to a quark in 
a higher Fock component (here Aqq) of the proton. 
Lower panel (c) : process where the photon cre- 
ates a qq pair leading to a non-diagonal transition 
between an initial 3g state and a final 4g q state. 



function to a 4qq state in the final wave function, as 
depicted in Fig. [8] (c). Such processes representing 
non-diagonal overlaps between initial and final wave 
functions are not positive definite and do not allow 
for a simple probability interpretation of the density 
p anymore. Only the processes shown in Figs. [8] (a) 
and (b), with the same initial and final wave function 
yield a positive definite particle density allowing for a 
probability interpretation. 

This relativistic dynamic effect of pair creation or 
annihilation fundamentally hampers the interpretation 
of density and any discussion of size and shape of 
a relativistic quantum system. Therefore, an inter- 
pretation in terms of the concept of a density re- 
quires suppressing the contributions shown in Fig. [8] 
(c). This is possible when viewing the hadron from 
a light front reference frame allowing for a descrip- 
tion of the hadron state by an infinite tower of 
light-front wave functions [21]. Consider the elec- 



tromagnetic (e.m.) transition from an initial hadron 
(with four-momentum p) to a final hadron (with 
four-momentum p') viewed from a light-front mov- 
ing towards the hadron. Equivalently, this corre- 
sponds to a frame where the hadrons have a large 
momentum-component along the z-axis chosen along 
the direction of the hadrons average momentum P = 
(p + p')/2. One then defines the light-front plus 
(+) component by a + = a? + a 3 , in a general four- 
vector a^, which is always a positive quantity for 
both quark or anti-quark four-momenta in the hadron. 
When we now view the hadron in a so-called Drell-Yan 
frame |22J, where the virtual photon four-momentum 
q satisfies q + = 0, energy-momentum conservation 
will forbid processes in which this virtual photon splits 
into a qq pair. Such a choice is possible for a space- 
like virtual photon, and its four-momentum or "virtu- 
ally" is then given by q 2 = —q 2 = — Q 2 < 0, where 
q± is the transverse photon momentum lying in the 
xy-plane. In such a frame, the virtual photon only 
couples to forward moving partons, i.e. only processes 
such as in Fig.[8](a) and (b) are allowed. We can then 
define a proper density operator through the + com- 
ponent of the four-current by J + = J° + J 3 [23] . For 
quarks it is given by 



<77 + <7 = 2 <2+<7 



+ • 



with 



-7 7+q, (5) 



where we introduced the q + fields through a field re- 
definition from the initial quark fields q involving the 
± components of the Dirac gamma matrices. The 
relativistic density operator J + , as defined in Eq. ([5]), 
is a positive definite quantity. For systems consisting 
of e.g. light u and d quarks, multiplying this current 
with the quark charges yields a quark charge den- 
sity operator given by J + (0) = +2/3 u(0)7 + u(0) — 
1/3 d(0)7 + (0)d(0). Using this charge density oper- 
ator, one can then define quark (transverse) charge 
densities in a hadron as |24l I25j: 



PoQ>) = 



x (P H 



d 2 q±_ 
(2n) 2 

' 2 ^ 



-i q ± -b 



1 

2P+ 



J+(0)|P+ - 



q± 



(6) 



where the hadron is in a state of definite (light-front) 
helicity. In the two-dimensional Fourier transform 
of Eq. (Jfij), the two-dimensional vector b denotes the 
quark position in the xy-p\ane relative to the position 
of the transverse centre-of-momentum of the hadron. 
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It represents the position variable conjugate to the 
hadron relative transverse momentum which equals 
just the photon momentum q±. 

The quantity po(b) has the interpretation of the 
two-dimensional unpolarized quark charge density at 
a distance b = \b\ from the origin of the transverse 
cm. system of the hadron. In the light-front frame, 
it corresponds to the projection of the charge density 
in the hadron along the line-of-sight. It is important 
to mind this difference to the interpretation in the 
non-relativistic case. 

The quark charge density in Eq. ^ does not fully 
describe the e.m. structure of the hadron, because 
we know that there are two independent e.m. FFs 
describing the structure of the nucleon. In general, 
a particle of spin S is described by (2S + 1) e.m. 
moments. In order to fully describe the relativistic 
structure of a hadron one needs to consider addition- 
ally the charge densities in a transversely polarized 
hadron state yielding a transverse charge distribution 
PTs±- We denote the transverse polarization direc- 
tion by S± = cos 4>s e x + sin e y . The transverse 
charge densities can then be defined through matrix 
elements of the density operator J + in eigenstates of 
transverse spin as [26l l27l 128] : 

PT (b)= f f^ e -i?x-bj_ 
PTsA°)-J (2 ^ )2 e 2p+ 

x(P + ,f ,sJJ+|P+,^, S± ), (7) 

where sj_ is the hadron spin projection along the di- 
rection of S±. Whereas the density p\ for a hadron in 
a state of definite helicity is circular symmetric for all 
spins, the density prs ± depends also on the orienta- 
tion of the position vector b, relative to the transverse 
spin vector S±. Therefore, it contains the information 
on the hadron shape, again projected on the plane 
perpendicular to the line-of-sight. 

The light-front wave functions and light-front den- 
sities discussed above are defined at equal light-front 
time (x + = 0) of their constituents. When con- 
stituents move non-relativistically, it does not make 
a difference whether they are observed at equal time 
(t = 0) or equal light-front time (x + — 0), since 
the constituents can only move a negligible small dis- 
tance during the small time interval that a light-ray 
needs to connect them. This is not the case, how- 
ever, for bound systems of relativistic constituents 
such as hadrons, or even when considering bound sys- 



tems of non-relativistic constituents if the reference 
system moves relativistically. Considering the latter, 
one may ask, for example, how the equal-time wave 
functions of bound state systems such as the hydrogen 
atom or positronium, which are non-relativistic bound 
states in QED at leading order in the fine structure 
constant a, transform when they move relativistically 
[29j [30] . The case of positronium has been studied in 
Ref. [29], by using the exact Bethe-Salpeter equation. 
It was found that the equal-time wave function of its 
leading |e + e~) Fock state contracts as expected from 
classical relativity. However, calculating the bound 
state in relativistic motion also necessitates the inclu- 
sion of an |e + e~7) component in the wave function, 
whose probability is of order 0(a). Unlike classical 
transformation laws, it was found that this photon 
amplitude, which can be considered as a quantum 
fluctuation, does not contract [29]. Therefore, the 
Fock components have a more complex shape change 
than just a flattening of their distributions along the 
direction of the boost. Only very recently, relativistic 
bound states, such as qq states in QCD have also been 
studied at lowest order in h [30J. It was found that 
the resulting equal-time wave functions have unique 
Lorentz transformation properties, ensuring the cor- 
rect dependence of the bound state energy on the 
center-of-mass momentum. A full understanding of 
the boost properties of bound state wave functions 
would allow to relate rest frame densities with light- 
front densities. For a relativistic many body system 
as a nucleon, for which this problem still needs to be 
fully solved theoretically, light-front densities provide 
at present a consistent theoretical method to image 
the charge density of quarks. Such imaging is anal- 
ogous to a flash photograph where different parts of 
the exposed object are reached by the light ray at 
different times. 

As summarized in the previous section, e.m. FFs of 
the nucleon are well measured experimental quanti- 
ties. We will, therefore, discuss the relativistic spa- 
tial shape as derived from these FFs. For a nu- 
cleon in a state of definite helicity, the transverse 
quark charge density is obtained from Eq. (^ by tak- 
ing the two-dimensional Fourier transform of its Dirac 
FF F 1 = (G e + tG m )/(1 + t) as [21 [25]: 

PQ (b) = j°° d ^QJ (bQ)F x (Q 2 ), (8) 

where J n denotes the cylindrical Bessel function of 
order n. Note that p only depends on b = \b\. 
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Figure 9: Quark transverse charge densities in the proton (left panels) and neutron (right panels). 
The upper panels show the density in the transverse plane for a nucleon with definite helicity. The 
lower panels for a nucleon polarized along the x-axis. The light (dark) regions correspond with largest 
(smallest) values of the density. For the proton e.m. FFs, the empirical parametrization of Arrington et 
al. [55] is used. For the neutron e.m. FFs, the empirical parametrization of Bradford et al. [33] is used. 



On the other hand, the information encoded in the 
Pauli FF F 2 = (G M - G E )/(1 + r) is connected 
to a nucleon in a transverse spin state. For a nu- 
cleon polarized along the positive ir-axis, the trans- 
verse spin state can be expressed in terms of the light 
front helicity spinor states by : \s± = +1/2) = 
(|A = +1/2) + e 4fe |A = -1/2)) /V2. The Fourier 
transform of the expression given in Eq. dTI in a state 
of transverse spin s± = +1/2 then yields [26j 



PTi,(b) 



sin 



PA (b) 

b - &s) 



dQ Q 2 
2tt 2M n 



Ji(bQ)F 2 (Q 2 ). 



(9) 



The second term, which describes the deviation from 
the circular symmetric unpolarized charge density, de- 
pends on the orientation of the transverse position 
vector b = b(cos<pbi x + sin^e^), relative to the 
transverse spin direction. 

In Fig. [9] the transverse charge densities in a nu- 
cleon, polarized transversely along the x-axis, i.e. for 
4>s = 0, are extracted based on the empirical informa- 
tion on the nucleon e.m. FFs which, however, does 
not yet contain the most recent data presented in 
Section[2] For the proton e.m. FFs, the empirical 
parametrization of Ref. [32J is used, whereas for the 
neutron e.m. FFs, the empirical parametrization of 
Ref. [33] is taken. One notices from Fig.[9]that polar- 
izing the proton along the x-axis leads to an induced 
electric dipole moment along the positive y-axis equiv- 
alent to the anomalous magnetic moment /j,^. This 
field pattern due to the induced electric dipole is a 
consequence of special relativity. The nucleon spin 
along the x-axis is the source of a magnetic dipole 
field B. An observer moving towards the nucleon 
with velocity v will see an electric dipole field pat- 
tern with E' = —j(v x B) giving rise to the observed 
asymmetry. 

For the neutron, one notices its charge density 
gets displaced significantly due to its large negative 
anomalous magnetic moment /ijv = —1.91 yielding 
an induced electric dipole moment along the negative 
y-axis. 

To see the impact of the most recent data, 
presented in Section[2] we compare in Fig. 10 the 
new data of Bernauer et al. [7] with the previous 
parametrization of world data, as performed by Ar- 
rington et al. [32J. In order to extract charge densities, 
one requires a form factor parametrization over all 



values of Q 2 . Because the Bernauer et al. data only 
provide a precision measurement of Ge p and Gmp 
for Q 2 < 0.4 GeV 2 , to fully quantify their impact on 
quark charge densities requires a new global analy- 
sis combining the previous data with these new data. 
Here we will perform a first estimate of this by using 
a parametrization which smoothly connects the new 
high precision data at low Q 2 and the Arrington et 
al. [32J parametrization at larger Q 2 . This interpo- 
lation function is displayed in Fig.[l0] and is used to 
extract the two-dimensional quark charge density in 
a proton in Fig.[TT] One readily sees that the new 
high precision data have a direct impact on the ex- 
tracted charge densities at large distances, typically 
larger than about 1.5 fm. By comparing the extracted 
density, using the previous fit to world data with the 
new fit, one sees that the new data lead to a sig- 
nificant reduction of the densities at distances larger 
than about 2fm. This is a direct consequence of the 
flatter behavior in Q 2 , for Q 2 < 0.3 GeV 2 , which the 
new data display for both Ge p and Gm p - 

In Fig. [12] we show the corresponding large distance 
behavior of the quark charge density in the neutron. 
The transition between the dashed blue curve and the 
solid red curve in Fig. [12] shows the impact of recent 
precision data at low Q 2 for the neutron FFs. These 
lead to a sizable enhancement in the extracted den- 
sities at distances larger than 1.5 fm. It is also of 
interest to compare these light-front densities with 
the static densities as discussed in section 2. For a 
non-relativistic system, one can extract from the 3- 
dimensional static density of Eq.Q, with intrinsic 
form factor p(k) = G*£;(fc 2 ), a 2-dimensional static 
density as : 



P2d(b) = 



+ 00 



dz p 3d (\/b 2 + z 2 ), 



+ 00 



dQ 
2tt 



QMbQ)G E (Q 2 )- (10) 



One notices that this static 2-dimensional density has 
the same form as the light-front density, see Eq. ([8]), 
with the crucial difference that in the static density 
the Sachs electric FF Ge appears, whereas for a rel- 
ativists system, the proper light-front charge density 
involves the Dirac FF F\. Since the large distance 
behavior is mostly impacted by the low Q 2 data, 
where Ge is dominated by F\, one expects a qual- 
itatively similar behavior at large distances between 



both pictures. This is illustrated in Figs. 11 and 12 
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Figure 10: Fits to the proton electric (upper panel) 
and magnetic (lower panel) FFs. Blue (dashed) 
curve: Arrington et al. [32 , black (dotted) curve: 
Bernauer et al. W] (only up to Q 2 = 0.4 GeV 2 ). The 
red curve presents a smooth connection between 
the Bernauer fit at low Q 2 and the Arrington et al. 
fit at larger Q 2 . 



where the light-front densities (solid red curves) are 
depicted along with the 2-dimensional static densi- 
ties (dotted black curves). One notices that for the 
proton, both densities approach each other at large 
distances pointing to a large tail in the charge dis- 
tribution. The corresponding picture for the neutron 
shows that both light-front and static densities dis- 
play a negative charge density for distances larger than 
about 1.6fm, which can be associated with a negative 
pion cloud in the outer region of the neutron. 

A combination of the FF data for the proton and 
neutron allows to perform a quark flavor separation 
and map out the spatial dependence of up and down 
quarks separately. The flavor separated FFs, invoking 



Figure 11: Large distance behavior of the unpolar- 
ized quark transverse charge density in the proton. 
The dashed blue curve uses the Arrington et al. 
parametrization [32) . The solid red curve shows the 
impact of recent high precision data at low Q 2 by 
using a smooth connection between the Bernauer et 
al. [7] fit at low Q 2 and the Arrington et al. fit at 
larger Q 2 (see solid red curves in Fig. 10 1. For com- 
parison, the solid green curve is the 2-dimensional 
projection of the static charge distribution accord- 



ing to Eq. (10), using the interpolating fit for Ge p 
(red curve in top panel of Fig. 10). 



isospin symmetry, are defined as 

Fl,2u = 2i*i .2 P + F12 n > 

Fl,2 d = Fl,2 p + 2^1,2 n - 



(11) 



For the Pauli FFs, it is convenient to divide out the 
normalizations at Q 2 = 0, given by the anomalous 
magnetic moments k u = 2n p + n n , Kd = k p + 2n n . 

We show the ratio of down/up flavor FFs in 
Fig. [13] and compare two empirical parametrizations 
with the result of a Generalized Parton Distribution 
parametrization of up and down quarks. For the em- 
pirical parametrizations, the dashed blue curves in 
Fig. 13 represent a previous fit to world data, us- 
ing the proton fit of Ref. [32J (dashed blue curve of 



Fig. 10 ) and the dipole type neutron fit of the Friedrich 



Walcher parametrization [16] (solid red curve in 
Fig.|6[). The solid red curve in Fig. 13 shows the im- 



pact of recent data at low Q 2 both for the proton [7j 
(solid red curve in Fig. 10), as well as the updated 
Friedrich-Walcher parametrization [16J for the neu- 
tron electric and magnetic FFs. One clearly sees from 



Fig. 13 that the down quark flavor FFs have a faster 
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Figure 12: Large distance behavior of the unpo- 
larized quark transverse charge density in the neu- 
tron. The dashed blue curve is the smooth part 
of the Friedrich-Walcher parametrization [TBi . The 
solid red curve is the updated Friedrich-Walcher 
parametrization which includes the recent neutron 
FF data at low Q 2 . For comparison, the solid green 
curve is the 2-dimensional projection of the static 



charge distribution according to Eq. ( 10 ), using the 
most recent fit for Geu- 



0.5 




2 

Q 2 [GeV 2 ] 



1.3 

,1.2 
< 

=1.1 



■o0.9 
X 

SO. 8 
0. 7 
0. 6 




Q 2 [GeV 2 ] 



fall-off than the up quark flavor ones. In the phe- 
nomenological GPD parametrizations this is encoded 
through a down-quark distribution which drops faster 
at large momentum fractions x than the up-quark dis- 
tribution. Using Eq. ([8]), we can then extract the ratio 
of up/down quark densities in the nucleon, which is 
displayed in Fig. [14] If the down and up quarks would 
have the same spatial distribution in the nucleon, the 
ratio as displayed in Fig. 14 would be one. We see 



however that in the center region of the proton, at 
distances smaller than about 0.5 fm, down quarks are 
less abundant than up quarks. The down quarks have 
a much wider distribution and are shifted to larger 
distances, dominating over up quarks between 0.5 to 
1.5 fm. At large distances, larger than about 1.5 fm, 
one clearly sees the impact of the recent data which 
results in a factor 2 change in the density as compared 
to previous fits to world data. Although the contribu- 
tion of the large distance region to the total charge 
is very small, the new data allow to precisely map 
out the charge densities in the region well beyond the 
confinement radius, where the charge density can in 
turn be interpreted as a measure of the contribution 
of the pion cloud. 



Figure 13: Ratio of down/up quark Dirac (upper 
panel) and Pauli (lower panel) FFs. For the Pauli 
FFs, the anomalous magnetic moments have been 
divided out. The dashed blue curves represent a 
previous fit to world data, whereas the solid red 
curves shows the impact of recent high precision 
data at low Q 2 as described in the text. The dotted 
black curves represent the results of a Generalized 
Parton Distribution parametrization [17| of up and 
down quarks. 

4 Conclusion 

We have presented two ideas about the long range 
structure of the nucleon. The first is nonrelativistic 
in terms of a "bare nucleon" plus a pion cloud, and 
the second relativistic in terms of quarks and gluons. 
One may be tempted to believe that the second is 
more fundamental since it uses the elementary fields 
of the standard model of particle physics. The quan- 
tum theory of quarks and gluons, QCD, describes a 
very large domain of strong interaction physics indeed. 
However, at sufficiently low energies, hadrons may 
be described by effective field theories formulated in 
terms of fields with discrete quantum numbers. These 
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Figure 14: Ratio of down over up quark densities in 
the proton according to the parametrizations and 
curve conventions of Fig. |13| 



fields may be viewed as elementary in a certain do- 
main of validity, i.e. sufficiently low energies here. 
One prominent example of such an effective field is 
just the pion in Chiral Perturbation Theory emerging 
as the Goldstone Boson of the spontaneous breaking 
of chiral symmetry of QCD. 

In fact, we are not able to devise a quantitative 
description of the nucleon-nucleon force in terms of 
quarks and gluons. On the other hand, the meson 
exchange idea allows for a very precise description of 
this force. Therefore, it may be futile to ask the ques- 
tion which description is more correct. Frequently in 
physics we have to be content with a model allowing 
a description in a limited domain and, following from 
this, limited predictive power. 

This sometimes confusing situation is also revealed 
by the two extreme reference frames in which we have 
considered the structure of the nucleon: the brick-wall 
system implying an infinitely heavy nucleon and the 
light-front frame implying a nucleon moving with ap- 
proximately the speed of light. As we demonstrated 
in both frames, the long distance structure of the nu- 
cleon reflects the physics of the pion cloud. With 
the recent experimental and theoretical advances we 
may have come closer to a common picture. However, 
whether we will ever be able to devise a "final theory" 
in terms of the elementary fields is an open question. 
Actually most physics is understood in terms of emer- 
gent effective degrees of freedom as the examples of 
condensed matter physics and nuclear physics show 
overwhelmingly. 
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